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ABSTRACT 



SN 2007od exhibits characteristics that have rarely been seen in a Type IIP 
supernova (SN). Optical V-hernd photometry reveals a very steep brightness de- 
chne between the plateau and nebular phases of ~4.5 mag, likely due to SN 
2007od containing a low mass of ^^Ni. The optical spectra show an evolution 
from normal Type IIP with broad Ha emission, to a complex, four component 
Ha emission profile exhibiting asymmetries caused by dust extinction after day 
232. This is similar to the spectral evolution of the Type Iln SN 1998S, although 
no early-time narrow (~200 km s"-*^) Ha component was present in SN 2007od. 
In both SNe, the intermediate-width Ha emission components are thought to 
arise in the interaction between the ejecta and its circumstellar medium (CSM). 
SN 2007od also shows a mid-infrared excess due to new dust. The evolution of 
the Ha profile and the presence of the mid-IR excess provide strong evidence 
that SN 2007od formed new dust before day 232. Late-time observations reveal 
a flattening of the visible lightcurve. This flattening is a strong indication of the 
presence of a light echo, which hkely accounts for much of the broad, underlying 
Ha component seen at late-times. We believe the multi-peaked Ha emission is 
consistent with the interaction of the ejecta with a circumstellar ring or torus 
(for the inner components at ±1500 km s""^) , and a single blob or cloud of cir- 
cumstellar material out of the plane of the CSM ring (for the outer component 
at -5000 km s~^). The most probable location for the formation of new dust is 
in the cool dense shell created by the interaction between the expanding ejecta 
and its CSM. Monte Carlo radiative transfer modehng of the dust emission from 
SN 2007od implies that up to ~4 x 10""^ Mq of new dust has formed. This is 
similar to the amounts of dust formed in other core-collapse supernovae such as 
SNe 1999em, 2004et, and 2006jc. 
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Subject headings: supernovae: individual (SN 2007od) — supernovae: general — 
circumstellar matter — dust, extinction 



-4- 



Introduction 



The discovery of massive amounts of dust in luminous infrared (IR) high-z galaxies 



seen less than 1 Gyr after the Big Bang (IBertoldi et al. 



20031 ). challenges our understanding 



of major dust contributors in the Universe. After only 1 Gyr, there has not been time for 
low-mass stars to form and evolve to the asymptotic giant branc h, which acc ounts for about 



19981 ). Therefore, it 



half of the dust production in the Galaxy at the current epoch (iDwekl 
has been suggested that core-collapse supernovae (CCSNe) may be significant contributors 
to the dust production in high-z galaxies. Surprisingly, dust mass estimates derived from 
the observed SEDs (spectral energy distributions) of recent SNe have been much lower than 
necessary to account for the dust observed in early-time galaxies. The sample of measured 
SNe is small but it appears that a typical Type II SN is producing only 10~^ — 10~^Mq 



of dust, well be 



early Universe (IKozasa et al. 



ow the 1 Mf;^ needed to account for t 



1989 



Todini fc Ferrara 



le lar ge amounts of dust seen in the 



have formed only a few 10 ^M© of dust (iKotak et al. 
that larger masses (1O~^M0) of dust can be derivec 



20091) 



2001 ). SN 2004et has been fo und to 



Sugerman et al. 



(120061) showed 



is considered to be clumpy and mostly silicate, but 



in the case o: 



Meikle et al. 



SN 2003gd if the dust 



( 120071 ) q uestioned these 



f indin gs and derive smaller dust masses (10 '^M(^ ) for S N 2003gd. Finally 



Ercolano et al. 



(120071 ). using similar methods to 



Sugerman et al. 



( 120061 ). did not find that large amounts of 



dust had condensed in the Si-poor envelope of SN 1987A. 

There are three distinct signatures of dust formation in CCSNe. These include a 
drop in visual brightness due to increased dust extinction accompanied by an infrared 
(IR) excess, and the appearance of asymmetric emission line profiles caused by the greater 
attenuation by dust in the ejecta of redshifted emission than blueshifted emission. These 



signatures usually 



( iLucy et al. 



1989 



appear a few 



Wooden et al. 



lundred days after the explosion and on l y SN 1987A 



19931 ). SN 2003gd (ISugerman et al. 



2006 



Meikle et al. 



20071 ). and SN 2004et fiSahu et all 120061 : iKotak et al.l 120091 ) have shown all three. 



A few CCSNe have shown evide nce for early dust fo rmation. SN 1998S formed 



dust sometime between days 140-268 (iLeqnard et al 



formation signatures between days 50-75 (IMattila et al 



2000h. while SN 2006ic sh owed dust 



2008 



Smith et al. 



20081 ). The SNe 



which show signs of early dust formation are all of type Iln, where the ejecta interacts with 
the circumstellar medium (CSM) of the S N. The dust may forni in a c ool dense shell (C DS) 



between the forward and reverse shocks (I Chevalier fc Fransson 



1994 



Pozzo et al 



20041 ) 



The formation of dust in the CDS allows for an additional avenue of dust formation not 
available in normal Type II SNe. The capability to form dust both within the SN ejecta 
and in the CDS suggests that these interacting systems have potential for forming dust in 
greater quantities than normal Type II SNe, although to date they have produced dust 
amounts comparable to normal Type II SNe. In this paper, we present new spectroscopic 
and photometric data covering 800 days of the evolution of SN 2007od, an unusual Type 
IIP SN that shows evidence for CSM interaction and early-time dust formation. 



Observations 



SN 2007od is located in UGC 12846 at a distance of 2 4.5 Mpc ( 



It was discovered on 2007 November 2 with R ~14.4 mag (IMikuz &: Maticic 



Immler &: Brown 



20071). 



20071 ). and two 



days later was confirmed spectroscopic ally as a normal Type I 



SN 1999em at 10 days past explosion (IBlondin fc Calkins 



SN most closely resembling 



20071 ). Based on this, we assume 



an explosion date of 2007 October 25 (JD 2454398) throughout this paper. The new 
observations reported in this paper are summarized in Table 1. We have obtained visible 
spectroscopy and photometry spanning days 57-811 as well as mid-infrared photometry on 
days 300, 455 and 667 of SN 2007od. 
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Optical imaging was obtained in Johnson-Cousins BVI with the SMARTS consortium 
1.3m telescope at Cerro Tololo Inter- American Observatory (CTIO), Chile. All images 
were pipeline reduced, shifted, and stacked. Imaging and spectra were also obtained with 
Gemini/GMOS North and South (GS-2008A-DD-2, GS-2008B-Q-45, GN-2008B-Q-245, 
GS-2009B-Q-1). The g'r'i' images were reduced and stacked using the IRAF gemini 
packa ge. The instrumen t al g^r^i^ magnitudes w ere transformed to standard Johnson- Cousins 



VRI flWelch et al. 



2007 



Fukugita et al. 



19961 ). For each night the transformation involved 



a least-squares fit with a floating zero point. Late time images were obtained with the 
Wide Field Camera (WFC) on HST/ACS using the F435W, F606W, and F814W filters. 
These images were pipeline reduced, including drizzling and cosmic ray re moval, and 



trans: 



Sirianni et al. 



'ormations to the Johnson- Cousins BVI system were accomplished using 
(j2005l ). The large bandpass of the F606W filter includes the 6563 A Ha emission, which 
is the main source of the late-time luminosity. We believe this causes the transformed 
Johnson V magnitude from HST to appear brighter than the V magnitude obtained from 
the Gemini transformations. In the Gemini filter set, r' contains the Ha flux, and is a 
very small component to the Johnson V transformation but a large component to the R 
transformation. Therefore the R magnitude from Gemini would contain the Ha emission 
whereas Ha would be present in the V magnitude from HST. 

A BVRI photometric sequence of tertiary standard stars (shown in Figure 1) was 
derived for the SN 2007od field. A table of magnitudes for these standards is located in 
Appendix A. The VRI lightcurves of SN 2007od are shown in Figure 2. Uncertainties for 
the Gemi ni images wer e calcu lated by adding in quadrature the transformation uncertainty 



quoted in 



Welch et al. 



(120071 ). uncertainty from photon statistics, and the zero point 



deviation of the standard stars f or each epoch . The HST uncertainties consist of the 



transformation uncertainty from 



Sirianni et al. 



( I2OO5I ) and the uncertainty from photon 



statistics added in quadrature. We used the narrow Na I D absorption seen in the day 
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10 spectrum of SN 2007od (described below) to estimate the foreground extinction. The 
equival ent width of the ble nded Na I D lines is 0.50 ±0.05A. This corresponds to A„ ~0.39 



±0.04 flBarbon et al 



19901 ). Since the value is small, the photometry presented in Figure 2 



and Table 2 have not been foreground corrected. 

For the first three Gemini/GMOS South epochs, three integrations of 900s were 
obtained in longslit mode using grating B600 and a slit width of 0'.'75. For the last 
two epochs on days 672 and 699, six integrations of 1800s each were obtained. Central 
wavelengths of 5950, 5970, and 5990 A for the first three GMOS-S epochs were chosen 
to prevent important spectral features from falling on chip gaps, and central wavelengths 
of 6350, 6370, 6390 A were chosen for the last two epochs to also include the [Ca II] 
AA7291,7324 emission. A 2x2 binning in the low gain setting was used. The one GMOS-N 
epoch used the R400 grating with central wavelengths of 5460, 5480, and 5500 A. All other 
settings were the same as in the early GMOS-S epochs. Spectra were reduced using the 
IRAF gemini package. The sky subtraction regions were determined by visual inspection, 
and included 30 pixels on either side of the SN. The spectra were extracted using 20 rows 
centered on the SN. The spectra from each individual night were averaged and have been 
corrected for the radial velocity of UGC 12846 (1734 km s^^). The spectra obtained on 
days 672 and 699 have been combined to increase the S/N. The spectra on days 10 and 50 
were taken with the FAST spectrograph operating on the 1.5-m Tillinghast telescope at the 
F. L. Whipple Observatory (Stephane Blondin, personal communication). The spectra are 
shown in Figures 3 and 4. 

We obtained three epochs of Spitzer/IRAC imaging and one epoch each with IRS 
PUI Blue and MIPS. See Table 1 and Table 3. The Spitzer IRAC (3.6, 4.5, 5.8, and 8.0 
/xm), IRS PUI Blue (16 /xm), and MIPS 24 fim images were mosaicked and resampled 
using standard MOPEX procedures to improve photometric quality. Aperture photometry 
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was performed utilizing the prescribed aperture corrections obtainable from the Spitzer 
Science Center website^ Figure 5 presents the SED of SN 2007od at three epochs, Spitzer 
(day 300)/Gemini (day 309), Spitzer (day 455-471) /Gemini (day 458), and warm Spitzer 
(day 667) /Gemini (day 672). The data have not been corrected for foreground extinction. 
Statistical uncertainties presented in the figure represent la errors. It should be noted that 
the SN was not detected in the MIPS and IRS/PUI images obtained on days 471 and 455, 
respectively. Consequently, these data are upper limits. Table 3 lists the measured mid-IR 
fluxes of SN 2007od. 



3. Discussion 



3.1. Lightcurve Evolution 



A typical Type IIP SN lightcurve shows a plateau for the first ~60-100 days, followed 
by a sharp drop in luminosity of ~1.5-3 mag as the SN transitions into the nebular 



phase, after which ; 



( IWoosley fc Weaver 



le lig htcurve exhibits an exponential decline powered by ^^Co decay 



19861 ). The lightcurve of SN 2007od (Figure 2) shows a normal plateau 
phase followed by an exceptionally large drop in brightness as it is entering the nebular 
phase. The last photometry before the gap was obtained on day 82 (Chornock 2009, 
personal communication). If we assume the end of the plateau phase took place shortly 
after that time and we extrapolate the slope of the ^^Co decay backwards into the gap 



(days 82-232) when no photometry is available, we find that t he V brightness 



minimum of 4.5 mag over the typical twenty d ay drop period ( 



mhamdi et al. 



'ades a 



2003ah . 



20091 ). have suggested 



Woosley fc Weaver! (119861 ) . and more recently (IKasen fc Woosleyi 
that if only a small amount of ^^Ni is synthesized in the SN explosion the lightcurve will 



^http: / / ssc.spitzer.caltech.edu/obs/ 
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show a steep drop between the plateau and nebular phases due to a lack of sufficient 
thermal energy. Generally, such a steep drop is only seen in low- luminosity (M^, ~-14) 
SNe such as 1999br, 1999eu and 2005cs, that have very low inferred ^^Ni masses (2 - 8 x 



lO'^ M^I flPastorello et al. 



122 ( ISollerman et al 



|2004J ). SN 1994W, which dropped ~3.5 mag from days 110 to 
19981 ) is the only other luminous (M^, ~-18) Type IIP besides SN 
2007od t o show such a steep dro p in brightness. Its ^^Ni mass was estimated to be ~2.6 x 



lO-^MfT 



0- 



SoUerman et al 



( 119981 ) suggested that after its 3.5 mag drop from the plateau, SN 
1994W declined much faster than the ^^Co decay once in the nebular phase, possibly due 
to dust formation. 

After day 458 the lightcurve of SN 2007od shows a deviation from ^^Co decay in the 
form of a flattening of the lightcurve between days 659 and 811. By day 811, assuming a 
constant brightness has been maintained from day 672, the difference in an extrapolated 
^^Co decay is then ~2.2 magnitudes in V. At late times, the Gemini and HST V magnitudes 
differ significantly due to differences in the filter bandpasses as discussed in Section 2. 
However, the constancy of the three epochs of HST/ACS V and / photometry from day 
659 to 811 show that the flattening of the lightcurve is real. It is likely being caused by 
a light pulse from the SN reflecting off of surrounding material and creating a light echo. 
Given the distance of the SN (24.5 Mpc), platescale of the WFC detector (0'.'05 pix^^) and 
FWHM of the stellar PSF (1.9 pix), a source could have a spatial extent on the plane of 
the sky of up to 0.12"or 46 It-yr and still be unresolved on the detector. In particular, we 
find through experimentation that a source with the flux of this SN and with a FWHM 
up to 1.2 times that of the PSF will not leave detectable residuals when subtracted away, 
therefore we allow t he physical size of the emitting region to be up to 0.12". Using the 



light echo equation ( ICouderc 



1939 



Sugerman 



20031 ). dust could be up to 475 It-yr in front 



of the SN on day 811 and still remain unresolved in the HST imaging. Thus, the echoes 
could be produced from either circumstellar or interstellar dust. Given that the foreground 
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extinction to the SN is low [Ay ~ 0.4), any dust producing a light echo should be optically 
thin. Typically, the contributi on from a light ec ho becomes significant when a SN has 



faded ~8 mag from maximum ( iPatat et al 



20061 ). SN 2007od was more than 10 mag below 



maximum light by day 672 when the contribution of the echo to the lightcurve became 
significant. A qualit ative c o mpar ison of the lightcurve and evolution at late times of SN 



2007od to models by 



Pat at 



(120051 ) suggest that the echo is more likely due to an interstellar 
sheet than circumstellar dust. We discuss possible contributions of an echo to the spectral 
evolution of SN 2007od in Section 3.2. 

As seen in Figure 2, the evolution of the V — I colors is normal for a dust forming 
Type IIP SN. We have only included the V — I colors from the Gemini observations due to 
the bandpass considerations of HST described in Section 2. At the time of dust formation, 
previous photometry of SNe 1999em and 2004et showed only small (<0.1 mag) changes in 
their V — I colors. The V — I color evolution of SN 2007od before and after the gap in 
observations (days 82-232) when we think dust formation began, is consistent with these 
other SNe. The lack of change in the V — I color evolution indicates that, in general, the 
color evolution is not sensitive to dust formation. This is likely caused by either a very 
small amount of dust being formed, or larger amounts forming in clumps. If the dust is 
clumpy, then the color evolution would not refiect the optical depth of dust present and the 
colo rs would not be a s red dened as they would be if the dust were uniformly distributed 



e.g. 



Witt fc Gordonl fll996[ ). 



3.2. Spectral Evolution 



The spectral evolution of SN 2007od shows behavior not previously seen in other 
SNe. Figure 3 shows the full optical spectrum of SN 2007od at 10, 232, and 458 days. 
Early-time spectra showed SN 2007od to be a normal Type II SN with a blue continuum. 
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broad Ha emission with a P-Cygni profile, and an expansion velocity of ~10000 km s ^ 



(IBlondin fc Calkins 



20071 Chornock 2009, personal communication). Other notable spectral 
features at early times include broad H/3 emission, weak Na I D AA5890,5896 absorption, 
and [O I] AA6300,6363 emission that was weaker than that normally seen in Type II SNe. 
There wa s no sign of either a narrow (FWHM ^ 200 krn s~^) Ha component nor X-ray 



emission (iBlondin fc Calkins 



2007 



Immler fc Brown 



20071 ). typically seen in a Type Iln SN. 



If present, the X-ray luminosity of SN 2007od was < 1.0 x 10^° erg s ^, sig nificantly smal ler 



than 



2007a 



uminosities measured in Type Iln SNe such as 2006jd and 2007pk (jimmler et al. 



By the time of the first Gemini/GMOS epoch (day 232) the spectrum of SN 2007od 
showed evidence for strong interaction with its GSM, an event usually associated with Type 
Iln SNe. The Ha emission profile had changed dramatically from the early-time spectra, 
with the single broad component evolving into at least three intermediate width (FWHM ~ 
1700 km s~^) emission components on top of a fainter, broader emission component (Figure 
4). Two of these intermediate width components are arranged symmetrically about zero 
velocity at ±1500 km s~^ and the third, which becomes more prominent at later epochs, lies 
at -5000 km s^^. We should note that the narrow emission at zero velocity is HII emission 
from the host galaxy not associated with the SN. Underlying the three intermediate width 
components, the broad Ha emission present in the early-time spectra can still be seen, 
and it is likely that it is responsible for the broad red wing of the Ha emission profile (see 
Figure 4). Broad, blende d Fe II emission between 5100-5400 A similar to that seen in SN 



1998S flPozzo et al. 



20041 ) appeared by day 232, and [Ga II] AA7291,7324 emission became 



prominent by day 458 (Figure 3). 

A similar mult i- component structure was seen in the Ha profile of SN 1998S. Unlike 
SN 2007od, SN 1998S was classified as a Type Iln, and showed both broad and narrow Ha 
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emission components at early-times due to the expanding ejecta and photo ionized CSM 



mate rial, respectively. The narrow component disappeared before day 108 flLeonard et al 



20001). and sometime before day 250, a triple-peaked Ha emission profile had developed 



( iGerardy et al. 



2000l ) (Figure 4). The three intermediate width components persisted in 



SN 1998S from approximately day 250 to day 660, exhibiting a progressive fading of the 



red and central components with respect to the 



was attributed to dust formation (IGerardy et al. 



jlue component; a phenomenon that 



2000l ). The timescale for the evolution 



of the SN 1998S spectrum is similar to that of SN 2007od, which showed a single broad 
Ha emission component from days 10-114, followed by a period (114-232 days) when the 
SN was unobservable, and finally the emergence of multiple-component emission by day 
232. In SN 2007od, the blueshifted component at -1500 km s~^ is much stronger than its 
corresponding component at -1-1500 km s~^. Like SN 1998S, the simplest explanation for 
this asymmetry is that newly formed dust is extinguishing the redshifted, far-side emission 
more than the blueshifted, near-side emission. The ratios of the inner peak heights would 
require an optical depth of approximately 1.0 at Ha for the asymmetry to be attributable 
to dust attenuation. If the dust existed in the CSM before the SN explosion, then it would 
lie outside of the ejecta and all of the emission components would be equally reddened. It 
is thought that the outer peaks in the Ha profil e of SN 1998S originated in the interactio n 



of the eiecta with a circumstel 



Pooley et al. 



2002 



Pozzo et al 



ar ring or torus flGerardv et al 



2004J : 



Fransson et al. 



200G 



Leonard et al. 



2000 



20051 ). If we assume the same is true 



of the asymmetric inner peaks of SN 2007od, then there must be dust between the source of 
the near-side (blueshifted) and the far-side (redshifted) emission components in the torus 



(Figur e 6). Moreover, if this dust is forming in the CDS, as was suggested by 



Pozzo et al. 



(120041 ) for SN 1998S, then the emission would likely be coming from a radiative forward 



shock (jPozzo et al. 



2004 



Fransson et al. 



20051 ) as one would expect equal attenuation of 



both redshifted and blueshifted emission originating in the reverse shock (see Figure 6). It 



- 13 - 



is also worthwhile to note that iKotak et al.l (|2009[ ) has recently alluded to CSM interaction 
producing a similar multi-component Ha profile in SN 2004et, but occuring at least 400 
days after it was seen in SNe 1998S and 2007od. 

In SN 2007od, the ±1500 km s~^ peaks appear sometime between day 114 and 232. 
Given the maximum velocity of the ejecta of 10000 km s~^, this constrains the distance 
of the inner radius of the CSM responsible for these profile peaks to ~700-1300 AU. 
Assuming a progenitor wind of 10 km s~^, this CSM was formed during a mass loss event 
ending 300-600 years prior to the SN explosion. With the emergence of the -5000 km s"^ 
components between days 114 and 309, a similar analysis places this CSM at 700-1700 AU 
and the formation period between 300-800 years prior to explosion. 

The appearance of the asymmetric components between days 114 and 232 implies that 
dust formation began much earlier than typically seen for non-interacting CCSNe. These 
SNe (SN 1987A, 2003gd, 2004et) showed developing as ymmetries between days 300 and 60 



post 



2006 



explosion due to dust condensation in the ejecta (ILucy et al. 



1989; 



Sahu et al 



Sugerman et al. 



20061). Furthermore, the [Ca II] A7291 emission profile in the day 458 
spectrum shows only the one broad component, with no intermediate-width components 
present (Figure 3). This implies that the complex emission seen in Ha is arising in the 
hydrog en-rich CSM lost fro m the star before the explosion and not in the metal-enriched 



ejecta ( IGerardy et al 



2000l ). These observation s are more consistent with dust formati on 



in the CDS as has been suggested in SN 1998S (jPozzo et al. 



2004 : 



Fransson et al. 



20051) 



The Ha profile of SN 1998S also showed a third peak near zero velocity that has been 

spherical distribut i on of shocked CSM 



clouds ( 


Gerarc 


V et al. 


2000: 


Fransson et al. 


2005 


). In the 



Leonard et al. 



2000 



Pooley et al. 



20021: 



Pozzo et al. 



2004; 



20051 ). In the case of SN 2007od, Figure 4 also shows the emergence of a 



third peak at -5000 km s ^ that was hidden by the broad ejecta emission on day 232. This 
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blue peak becomes increasingly more prominent as the broad underlying ejecta emission 
begins to fade over time. By day 699, when we believe the light echo is contributing 
significantly to the flux, this feature is still prominent. Because a light echo depicts only 
the early-time emission, it is not possible to explain the -5000 km s~^ peak with a light 
echo scenario. The nature of this third peak represents the starkest difference between SNe 
2007od and 1998S. As previously mentioned, in SN 1998S the peak not associated with 
the CSM torus was observed at zero velocity and attributed to a spherical distribution of 
shocked CSM clouds. Although a similar solution can not be invoked to explain a peak 
centered at -5000 km s~^ as seen in SN 2007od, it is possible that this peak comes from the 
emission of a single cloud or blob of shocked CSM. If this scenario is correct, because the 
torus and the blob are seen at different velocities, the blob cannot be located in the plane 
of the torus. At the same time, the blob would be required to have a velocity component 
along our line of sight larger than that of the torus. Therefore, if the torus is nearly edge on 
as viewed from Earth (Figure 6), as hypothesized in SN 1998S, the blob cannot be moving 
orthogonal to the plane of the torus as this would severely limit its velocity component 
along our line of sight. On the other hand, if the inclination^ of the ring is less than 90°, 
the blob could be moving orthogonally, but it would require a sufficiently thick torus to 
retain the dust extinction necessary to explain the asymmetric inner peaks. The scenario 
presented here to account for the observations is closely based on the scenario used to 
explain the observations of SN 1998S. This torus -|- blob scenario is introduced in order 
to illuminate the manner in which the multi-component Ha emission may arise, but the 
existing data do not allow the detailed geometry of SN 2007od to be ascertained. 

From days 232 to 699, the relative strength of the emission component at -5000 km s~^ 
is increasing as the inner pair of emission components (± 1500 km s~^) slowly disappear. 



2 2 — 90° is edge-on 
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Again if we assume the SN 1998S scenario, then this is explained by the ejecta moving 
through, and completely destroying, the torus but not the blob. The underlying broad 
component is also present at all epochs. The broad component is likely a combination 
of the residual ejecta emission and, given the flattening of our lightcurve ~10 mag below 
maximum, a light echo. The contribution to the total broad component by the light echo 
is expected to grow with time as the ejecta emission fades. Based on the lightcurve, on 
day 238 the light echo flux would be ~ 20x fainter than the directly transmitted ejecta 
emission, and on day 348 the echo would still be ~8x fainter. However, by day 672/699, 
the echo flux will likely be greater than the directly transmitted ejecta emission. Therefore 
the echo should not be a major contributor to the Ha proflle at early times, but would be 
a major contributor of the late-time proflle. Figure 4 shows the flts to the day 232 and 
348 Ha proflles. We have used Gaussians to approximate the intermediate width peaks 
and a scaled day 50 spectrum to approximate the residual ejecta emission. We should note 
that the scaled day 50 spectrum is only an approximation to the day 232 and 358 ejecta 
emission. Realistically, in addition to fading over time, the proflle would likely lose its 
P Cygni absorption, and its overall width would likely decline. The flat-topped proflle of 
the -5000 km s~^ component is shown to be well flt by two intermediate width Gaussians 
possibly due to density enhancements within the GSM cloud or blob. At later times the 
proflle becomes too complex to flt. On days 458 and 672/698, the -5000 km s~^ peak 
begins to dominate and its two component nature becomes more clear. Moreover, there is 
a continuing fading of the inner peaks to the point where they disappear into the broad 
emission by day 672/698. As previously mentioned, by this last epoch, the broad emission 
is a combination of the light echo spectrum and the weak residual ejecta emission. However, 
the residual flux beyond 5000 km s~^ is unlikely to originate from either of these two sources 
and remains unexplained. There is also an indication that the radial velocities of both sets 
of components get smaller with time. This is similar to the behavior of the blue component 
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of SN 1998S which moved to smaller velocities from day 249 to da y 658, and was exp lained 



as an interaction of the shock with more massive sections of CSM (jPozzo et al. 



2004h . 



3.3. Radiative Transfer Modeling 



The SEDs of SN 2007od shown in Figure 5 all show strong IR excesses attributable 
to warm dust. The SED at 300 days can be fitted well by the sum of a blackbody with a 
temperature of 5100 K representing hot, optically-thick gas in the ejecta and a modified 
blackbody with a temperature of 580 K subject to a emissivity law representing the 
dust. A similar fit to the day 455 SED uses a 5100 K blackbody and a 490 K modified 
blackbody, and for day 667 a 5100 K blackbody and a 600 K modified blackbody. The 
fits are shown in Figure 5 and the relevant parameters are compiled in Table 4. The IR 
emission could arise from pre-existing dust in the CSM, and/or new dust condensing in 
either or both of the CDS and the ejecta itself. IR echoes from pre-existing CSM dust 
heated by the early-time luminosity are often seen in other SNe. The position of the dust 
in SN 2007od cannot be inferred from IR observations alone, but when combined with the 
evidence provided by the asymmetry in the Ha emission components, it suggests that the 
IR emission is coming from dust forming in the ejecta or the CDS, and not from an IR 
echo. The very early onset of dust formation (between 114-232 days) further suggests that 
the dust formation is occurri ng in the CDS sinc e the ejecta at this time would still be too 



warm for dust condensation (IKozasa et al. 



199lh . 



In order to better understand the amount and composition of the dust forming 
within SN 2007od we used our 3D Monte Carlo radiative transfer (RT) code MOCASSIN 



(lErcolano et al. 



2005 



and references therein). For the purposes of our RT models we will 
assume the dust and source luminosity is spread throughout a spherical, expanding shell 
with an inner radius, Rj„, and an outer radius radius. Rout- The luminosity at each location 
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is proportional to the local density. With this simple model, we can quantify the total 
amount of dust in the system as a whole (ejecta + CDS) and monitor its variation from day 
300 to 667. We have modeled two different dust distributions. Following radiative transfer 



models of other Type II SNe (ISugerman et al. 



2006 



Ercolano et al. 



20071 ). the "smooth" 



model for SN 2007od has the dust uniformly distributed throughout the shell according 
to a r~^ density profile, while the "clumpy" model reflects an inhomogeneous distribution 
of spherical clumps embedded within an interclump medium of density p. The irradiating 
photons are only produced in t he interclump medium, while the clumps are assumed to be 



dark (see 



Ercolano et al. 



( 120071 ) for further discussion of this issue). We have adopted a 



standard grain siz e distribution of a between 0.005 and 0.05 /im (ISugerman et al. 



Meikle et al. 



2007 



Kotak et al. 



2006; 



20091 ) and have tested two different concentrations, i.e. a 
dust model dominated by silicates (75% Si; 25% AC) and one dominated by amorphous 
carbon (75% AC; 2 5% Si). For each of t hese models, we used the optical grain constants of 



Hannerl (119881 ) and 



Draine &: Led (Il984j ). respectively. 



The main input parameters for the model are the ejecta temperature and luminosity, 
the inner and outer radii of the shell, and the mass of dust present. Initial values of the 
luminosity, temperature, and approximate size of the dust shell are based on the results of 
our blackbody analysis above. Tuning the dust mass, along with slight variations in the 
other parameters, allowed us to obtain a good fit to the optical and Spitzer/IRAC points 
for the first two epochs. We should point out that without the 5.8 and 8.0 /xm data on day 
667, we were unable to make reliable constraints on the dust mass at this epoch. The best 
fit to the data was obtained with a smooth dust distribution and is shown in Figure 5 and 
the corresponding input parameters are listed in Table 5. This table also lists the dust mass 
calculated by our clumpy model. 



Although no attempts were made to fine tune the models to match the 16 and 24 ixm 
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flux upper limits on day ~455, care was taken to ensure that our preferred model did not 
predict emitted fluxes that would have been detected given our experimental setup. We 
found that we were unable to account for this relatively low flux at 16 and 24 /im with 
compositions dominated by silicate dust (see Figure 5). All other parameters being equal, 
transitioning from amorphous carbon dust to silicate dust has the effect of broadening 
the IR peak as the SED will then be increased by strong silicate emission features at 10 
and 18 /xm. These two emission bumps are clearly seen in the Si model in Figure 5. The 
conclusion is that the silicate dominated model predicts too much flux beyond 10 /xm at 
day 455 leading us to favor the amorphous-carbon dominated model with 75% AC and 25% 
Si. Given this concentration, our best smooth model predicts a total of 1.7 x 10^^ Mq and 
1.9 X 10~^ Mq of dust on days 300 and 455, respectively, with an increase in from 1.5 to 
2.0 over the two epochs. This implies a total of 2 x 10^^ Mq of new dust forming between 
day 300 and 455 after explosion. The values of r at 6600 A predicted by our MOCASSIN 
fits were 1.3 and 1.8 for days 300 and 460, respectively. These opacities are adequate to 
explain the inner peak height ratios discussed in Section 3.2. Although the constraints were 
less reliable on the dust mass at day 667, we were able to show that the available data is 
consistent with our previous epochs with a predicted dust mass of 1.8 x 10""^ Mq and 
= 2.75. The best clumpy model gave a total dust mass of 2.5 x 10~^ Mq and 4.2 x 10~^ 
Mq on days 300 and 455/667, respectively, implying the formation of 1.7 x 10~^ Mq of 
new dust formed. In each case, the quantity of dust formed is several orders of magnitude 
l ower than that required to contribute significantly to the dust budget of the early universe 



(IMorgan fc Edmunds 



20031). 
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4. Summciry 

SN 2007od is a very unusual Type IIP SN that shows strong evidence for early-time 
(<232 days) dust formation in the cool dense shell (CDS) created by the interaction 
between the ejecta and the circumstellar medium (CSM). Observations prior to day 82 for 
the optical lightcurve and prior to day 114 for the spectral evolution indicate SN 2007od 
was a normal Type IIP SN. Photometry and spectra obtained on day 232 revealed that a 
large post-plateau luminosity drop (~ 4.5 mag) and the development of multi-peaked Ha 
emission had occurred during the period the SN was unobservable. The sharp luminosity 
decline between days 82 and 232 may have been caused by a low mass (<3 x 10~^Mq) of 
^^Ni in the ejecta combined with the formation of new dust in the CDS. The presence of 
dust is supported by the red/blue asymmetries seen in the Ha; components as well as the 
presence of a large IR excess. In SN 2007od, the dust formed between days 114-232 which 
is much earher than seen for other Type IIP SNe, but in line with SNe that are interacting 
with the CSM such as SN 1998S. There is also an indication in the observed flattening 
of the late-time lightcurve that a light echo from surrounding circumstellar or interstellar 
medium is present around SN 2007od. Flux from this echo also would make a significant 
contribution to the Ha emission profile at late times. 

SN 2007od is unusual in that, while it is classified as a Type IIP, there is evidence for 
CSM lying close to the progenitor, even though it showed neither a narrow Ha component 
nor X-ray emission which are both present in typical Type Iln SNe. We believe the 
multi-peaked Ha emission in SN 2007od arises from a similar mechanism to that seen in SN 
1998S, and is consistent with the interaction of the ejecta with a circumstellar torus and 
blob. Assuming that scenario, the ± 1500 km s^^ emission arises from the ejecta interacting 
with a circumstellar torus, and the -5000 km s~^ arises from a single cloud or blob of 
shocked CSM out of the plane of the CSM. The use of a torus and blob is to illustrate the 
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ways that the multicomponent Ha emission may arise but should not be construed as the 
"true" geometry of the SN 2007od which is not well constrained by the data. 

Our radiative transfer models have indicated that up to 4.2 x 10~^Mq of dust has 
formed in the CDS of SN 2007od. The models indicate that 2.5 x lO^^M© was present 
on day 300 and an additional 1.7 x 10~^Mq formed by day 455. This amount of new 
dust is similar to estimates for other C CSNe. For example SN 1987A produced 7.5 x 



^ M(7) of dust (Ercolano et al. 



(IMattila et al. 



20071 ) and SN 2006jc an estimated 3 x lO^^M© of dust 
20081 ). Dust est imates for SN 2003gd currently yield values b etween 2 x 



lO-^M© and 2 x lO'^Mf;, from 







Sugerman et al. 



(120061 ) and 4 x lO'^M© from flMeikle et al 



20071 ). both based on the same data. There seems to be a significant variation in the 
amount of dust produced from SN to SN, but none of them approach the ~1 Mq of 
dust per SN needed to ac count for the mass i ve amounts of dust see n in early-time 



galaxies at high red-shift (IKozasa et al 



1989 



Todini fc Ferrara 



200ll ). Therefore, in the 



few dust-forming Type II SNe that have been observed, there is no evidence that these 
objects are major contributors to the dust-rich early galaxies seen at high redshift. The 
sample is still very meager, and more such SNe need to be observed with good tem- 
poral and wavelength coverage. Monitoring of SN 2007od is continuing in the optical and IR. 
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Table 1: Observation Summary 



Day 


JD 


Telescope 


Instrument 


Exposures 


Time/Exposure 


57 


2454455 


SMARTS 1.3m 


ANDICAM 


5 


20 


232 


2454630 


Gemini South 


GMOS Spectra 


3 


900 


232 


2454630 


Gemini South 


GMOS Imaging 


3 


20 


300 


2454698 


Spitzer 


IRAC 


12 


100 


309 


2454707 


Gemini South 


GMOS Spectra 


3 


900 


309 


2454707 


Gemini South 


GMOS Imaging 


1 


60 


348 


2454746 


Gemini South 


GMOS Spectra 


3 


900 


348 


2454746 


Gemini South 


GMOS Imaging 


2 


60 


455 


2454853 


Spitzer 


IRS PUI 


9 


15 


458 


2454856 


Gemini North 


GMOS Spectra 


3 


900 


458 


2454856 


Gemini North 


GMOS Imaging 


1 


60 


461 


2454859 


Spitzer 


IRAC 


12 


100 


471 


2454869 


Spitzer 


MIPS 


16 


30 


659 


2455057 


HST 


ACS/WFC 


4 


98 


667 


2454865 


Spitzer 


IRAC 


12 


100 


672 


2455070 


Gemini South 


GMOS Spectra 


6 


1800 


672 


2455070 


Gemini South 


GMOS Imaging 


2 


60 


699 


2455097 


Gemini South 


GMOS Spectra 


6 


1800 


709 


2455107 


HST 


ACS/WFC 


4 


98 


811 


2455109 


HST 


ACS/WFC 


4 


98 



(s) 
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Table 2: Optical Photometry of SN 2007od 



Day 


V 


R 


/ 


57 


14.86 ± 0.04 




14.17 ± 0.03 


232 


21.32 ± 0.04 


20.43 ± 0.02 


20.21 ± 0.03 


309 


22.05 ± 0.05 


21.36 ± 0.03 


21.10 ± 0.05 


348 


22.44 ± 0.05 


21.80 ± 0.03 


21.54 ± 0.05 


458 


23.51 ± 0.12 


22.66 ± 0.07 


22.77 ± 0.15 


659 


23.64 ± 0.18 




24.08 ± 0.35 


672 


24.67 ± 0.15 


23.82 ± 0.10 


24.10 ± 0.25 


709 


23.68 ± 0.21 




24.22 ± 0.50 


811 


23.76 ± 0.20 




24.21 ± 0.44 



Note. — Day 57 is SMARTS/ANDICAM photometry, and days 659, 709, and 811 are HST/ACS. The 
remaining epochs are Gemini/GMOS photometry. 
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Table 3: Spitzer Photometry of SN 2007od 



Day 3.6/um 4.5//m 5.8/um 8.0/um 16/um 24//m 

300 110.1 ± 1.6 156.4 ± 1.8 200.0 ± 7.6 199.5 ± 12.2 

461 37.9 ± 1.4 63.4 ± 1.6 90.0 ± 7.5 122.2 ± 11.8 <59 <35 

667 23.5 ± 0.7 33.3 ± 0.8 - - - - 



Note. — Fluxes are reported in /uJy. The 16/im photometry is from day 455 and the 24/im photometry 
from day 471. The day 667 IRAC photometry was obtained with warm Spitzer, which only operates in 3.6 
and 4.5 /xm. 

Table 4: Blackbody Fits 

Blackbody Optical Mod. Blackbody IR 

Epoch T (K) R (AU) L (Lq) T (K) R (AU) L (Lq) Uot. (Lq) 

300 d 5100 5.7 9.2e5 580 1824 2.4e6 3.3e6 

455 d 5100 2.9 2.3e5 490 1860 l.le6 1.3e6 
667 d 5100 1.6 7.5e4 600 767 5.1e5 5.8e5 



Table 5: Monte Carlo Radiative Transfer Models 



Smooth Clumpy 



Epoch 


AC/Si 


T(K) 


Rm (AU) 




Ltot. (L©) 




Md (Mo) 


Md (Mo) 


300 d 


0.75/0.25 


6000 


50.1 


2005 


3.4e6 


1.50 


1.7e-4 


2.5e-4 


455 d 


0.75/0.25 


6000 


50.1 


1537 


1.4e6 


2.00 


1.9e-4 


4.2e-4 


667 d 


0.75/0.25 


6000 


50.1 


936 


8.1e5 


2.75 


1.8e-4 


4.2e-4 
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Fig. 1. — The field of SN 2007od showing the location of the tertiary standards (designated 
alphabetically) listed in Table 6. The location of SN 2007od is marked. North is up and 
east is left. The image is 4.7' x 3.1' and was taken with Gemini/GMOS-S on 2008 August 
29 (day 309) with the g' filter. 
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Fig. 2. — Left: Johnson- Cousins VRI lightcurves of SN 2007od. R points are shown 
as circles and have been shifted up 2 mag, V are shown as squares, and I as triangles 
shifted down 3 mag. The filled symbols are from Gemini/GMOS and the open symbols from 
HST ACS/WFC. Additional photometry of SN 2007od obtained with no filter is plotted 
(asterisks). The uncertainties are listed in Table 2, and are in most cases smaller than 
symbol size. The dashed line represents the ^^Co decay. Note the drop of ~4.5 mag in V 
in 20 days, which is large for a Type IIP SN. The fiattening of the lightcurve after day 600 
shows the presence of a light echo which is ~10 mag fainter than the SN at maximum.The 
disagreement in V between Gemini and HST at late epochs is most likely due to filter 
bandpass differences (see text). Right: Evolution of Gemini V — I of SNe 2007od (filled 



squar es), 1999em (triangles), and 2004et (diamonds) (lElmhamdi et al. 



2003b 



Sahu et al. 



20061 ). The V — I colors of SNe 1999em and 2004et have been corrected for foreground 



extinction. 
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Fig. 3. — Optical spectra of SN 2007od on days 10, 232 and 458. All spectra have been 
normalized to the continuum and are shifted vertically by a constant. See text. 
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Fig. 4. — Left: Spectral evolution of the Ha profile of SN 2007od. The -5000 km s ^ 
component is easily visible by day 309, and the inner ±1500 km components disappear 
by day 672/6 99. 



(IPozzo et al. 



or comparison, the Ha profile of SN 1998S (day 305) has been included 
20041 ). Right: Deconvolution of the Ha profiles on day 232 and 309 are shown 
on the right side panels, and the fit is shown as a black line over the spectrum on the left side 
panels. The dashed-dotted line is a scaled day 50 profile, showing that the broad underlying 
emission at these epochs is consistent with being produced by directly transmitted light from 
the expanding ejecta. The intermediate width ± 1500 km s~^ and -5000 km s~^ components 
are typically produced by CSM interaction as shown by comparison with the SN 1998S 
spectrum. The double gaussians required to fit the -5000 km s~^ component likely represent 
density enhancements in the CSM. The day 672/699 spectra, shown on the left, is likely to 
be mostly due to a light echo which only became a significant contributor when the SN had 



faded by ~ 10 mag. 
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Fig. 5. — SED analysis of SN 2007od on days ~300 (diamonds), ~460 (triangles), and 
~667 (squares). Visible wavelengths were obtained with Gemini/GMOS and IR wavelengths 
with Spitzer/IRAC, IRS PUI, and MIPS. Dust emission is clearly visible in the 3.6-8 iim 
wavelength range. Left: SEDs with best fit blackbody curves. The blackbody data is given 
in Table 4. Optical data was fitted with a simple blackbody while the IR data was fitted 
with a modified blackbody subject to a emissivity law. Right: The corresponding best 
fit MOCASSIN dust models for each epoch. The best fits refiected a smooth distribution 
of 75% AC dust (solid lines). For comparison, the best fit 75% silicate model for day 455 
(dashed line) is shown and predicts too much fiux past lO/itm. MOCASSIN input parameters 
are shown in full in Table 5 for each model. 
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Fig. 6. — A ca r toon s howing; the geometry of SN 2007od (adapted from 



Mattila et al. 



and 



Fox et al 



(120081 ) 



( 120091 ) )■ The inner ±1500 km s^^ components may come from the forward 
shock created when the SN ejecta plowed into a central CSM ring, and the outer -5000 km 
s~^ component from the ejecta interacting with a blob of CSM out of the plane of the central 
ring. If the dust is forming in the CDS between the shocks, then the emission coming from 
the forward shock at the CSM ring explains why the near-side emission is less reddened than 



the far-side emission. See text. 
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A. Tertiary Standards for UGC 12846 

The photometric sequence of tertiary standard stars used for both SMARTS and 
Gemini photometry (Table 6) was obtained using the Y4KCam CCD on the Yale 1.0m 
telescope operated at the Cerro Tololo Inter-American Observatory on the night of 2008 
October 28. The YlKCam offered a field of view of ~20'x20' and employed a set of 
standard Johnson-Kron-Cousins broadband BVRI filters. The processing of all images 
to remove instrumental signatures was accomplished using the standard techniques of 
subtracting a median filtered bias frame and dividing by master twilight fiat field images 
for each filter. Also, the fringing effects in the I-band exposures were removed by scaling 
and subtracting a master fringe frame from the program images. In order to transform 



of the instrumental magnitudes to the standard sys tem, we o 



Dserved between 30 and 35 



different stars (depending on filter) from the lists of iLandoltl (120091 ) in addition to the SN 
2007od field. Importantly, these stars were selected due to their ample range in color, and 
they were observed both near zenith and at a high airmass (~1.75). Thus, they provided 
a viable sample of stars from which to derive accurate color and extinction coeffecients as 
well as magnitude zero-points for each filter. The resulting RMS scatter of the residuals 
from the transformations revealed that the photometry for selected stars in the SN 2007od 
field derived from this night's observations is accurate on the 1-2% level depending on filter. 

Table 6: Tertiary BVRI Standards for UGC 12846 



Star 


B 


V 


R 


I 


A 


16.989 ± 0.014 


15.938 ± 0.011 


15.338 ± 0.006 


14.820 ± 0.017 


B 


17.480 ± 0.020 


16.357 ± 0.011 


15.720 ± 0.006 


15.153 ± 0.018 


C 


17.266 ± 0.015 


16.109 ± 0.011 


15.403 ± 0.006 


14.786 ± 0.017 


D 


17.966 ± 0.019 


17.133 ± 0.014 


16.672 ± 0.010 


16.343 ± 0.068 


E 


18.146 ± 0.021 


17.240 ± 0.014 


16.708 ± 0.010 


16.267 ± 0.039 



